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CM We use QCD sum rules to study the recently observed meson Z+(4430), considered as a D* Di 

molecule with J p = 0~. We consider the contributions of condensates up to dimension eight and 
work at leading order in a s . We get mz = (4.40 ± 0.10) GeV in a very good agreement with the 
experimental value. We also make predictions for the analogous mesons Z a and Zbb considered as 
DtD\ and B* B\ molecules respectively. For Z s we predict mz s = (4.70±0.06) GeV, which is above 
the D*D\ threshold, indicating that it is probably a very broad state and, therefore, difficult to 
observe experimentally. For Zbb we predict mz bb = (10.74 ± 0.12) GeV, in agreement with quark 
model predictions. 
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In the last years many new mesons have been observed by the BaBar, BELLE, CLEO, DO and FOCUS 
collaborations. Among these new mesons, some have been considered as good candidates for tetraquark 
states like the D sJ (2317) Q, the D sJ (2460) 0, the X(3872) % and more recently the Z+(4430) Q. 
While there are many indications that the charmed mesons, D s j(2317) and D s j(2460), are not four-quark 
states [1101, this is not the case for the charmonium like states, A"(3872) and Z+(4430). The A(3872), 
| with quantum numbers J = 1 ++ , does not fit in the charmonium spectrum and presents a strong 

isospin violating decay, disfavoring a cc assignment. The Z + (4430), recently observed in the Z + — » iJ/t: + 
decay mode [3| , is the most interesting one since, being a charged state, it can not be a pure cc state. 

There are already many theoretical interpretations for the structure of the Z + (4430) meson: molecular 
D*Di state Q, tetraquark state [§, 9j, or a cusp in the D*D 1 channel [l(|. In ref. [g|, the authors have 
F-- interpreted the Z+(4430) meson as the first radial excitation of the diquark-antidiquark [cu][cc?] state, 

with J PC = 1^ . The low lying tetraquark state, [cw,][cd], is interpreted as the charged partner of the 
X(3872) meson. Supposing that the mass of the low lying tetraquark state, [cw][cc?] is compatible with 
the mass of the X(3872) meson, the mass difference between the Z + (4430) and the A"(3872) would be 
close to the mass difference between the ip' and J/ip: m^* — mjw = 590 MeV. Therefore, they arrive at 
mz ~ 3872 + 590 ~ 4460 MeV, which is compatible with the observed mass. In ref. Q, the closeness of 
the Z+(4430) mass to the threshold of D* + (2010)51 (2420) lead the authors to consider the Z+(4430) as 
a D*Di molecule. In this case, the allowed J p of Z would be 0~, 1~ or 2 _ , although the 2~ assignment 
is probably suppressed in the B — > Z (4430) K decay, by the small phase space. Among the remaining 
possible 0~ and 1~ states, the former will be more stable as the later can also decay to DD\ in s-wave. 
Moreover, one expects a bigger mass for a J p = 1~ state as compared to a J p — _ state. Therefore, 
in this work we use QCD sum rules (QCDSR) [ll], [l2|, [l^|, to study the two-point function of the state 
Z+(4430) considered as a D*Di molecule with J p = 0~. 

In a previous calculation, the QCDSR approach was used to study the X(3872) meson, considered as 
a diquark-antidiquark state, and a good agreement with the experimental mass was obtained 15]. If we 
suppose, as in ref. [8j, that the Z + (4430) is related to the first radial excitation of the X(3872), in the 
QCDSR approach its mass would be given by y/so, where sq is the continuum threshold. In ref. [6j it 
was found that a/so = (4.3 ± 0.1) GeV also in a very good agreement with the experimental mass of 
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Z + (4430). However, this is not a precise determination of the mass of the first excited state, since the 
continuum threshold gives only a lower bound for the mass of the first excited states. 

Considering Z + (4430) as a D*D% molecule with J p = CP, a possible current describing such state is 
given by: 

3 = -J= [(dal^icb^^Ub) + (d a 7 /1 75C )(C67' i Ub)] , (1) 

where a and b are color indices, We have considered the symmetrical state D* + + D*°D^ because it 
has positive G-parity, which is consistent with the observed decay Z + (4430) — > ■0 / 7r + . 
The two-point correlation function is given by: 

n(g)=t / d A x e**Q\T\j{x)3\tim- ( 2 ) 



On the OPE side, we work at leading order in a s and consider the contributions of condensates up to 
dimension eight. We calculate the light quark part of the correlation function in the coordinate-space and 
we use the momentum-space expression for the charm quark propagator. The resulting light-quark part 
is then Fourier transformed to the momentum space in D dimensions and it is dimensionally regularized 
at D = 4. 

The correlation function in the OPE side can be written as a dispersion relation: 



roc OPE(„\ 

n OPi V)=/ da? \p. + n mix ^(q 2 ) , (3) 

J 4m 2 s-q 2 



where p OPE (s) is given by the imaginary part of the correlation function: np OPE {s) — lm[Il (s)]. We 
get: 



p OPE (s) = pP ert (s) + p<«>( 8 ) + p^\s) + p mix (s) + p^\s) , (4) 



ith 



0n 



p<™>(«) = o, 



l-a 



p(G ^ {s)= rn^G)_ da , ^ _ ^ _ ^ [{a + &)< _ 



2 8 7T 6 

ct m in 0mi„ 

p mix (s) = 0, 

U mix{qq) ^ = m 2 c (qga.Gq) (qq) f 1 ^_ a(l - a) 



2 3 7r 2 J to 2 — a(l — a)q 2 



m 2 . — a(l — a)q 2 



(5) 



where the integration limits are given by a m i n = (1 — yl — 4m 2 /s)/2, a max = (1 + y/l — 4to 2 /s)/2 and 
Pmin = am 2 /(sa — m 2 ). The contribution of dimension-six condensates (g 3 G 3 ) is neglected, since it is 
suppressed by the loop factor l/16n 2 . In Eq. ^ the Yl mlx (vi) (g 2 ) term is treated separately because its 
imaginary part is proportional to delta functions and, therefore, can be easily integrated 

It is very interesting to notice that the current in Eq. ([T]) does not get contribution from the quark 
and mixed condensates. This is very different from the OPE behavior obtained to the scalar-diquark 
axial-antidiquark current used for the X(3872) meson in ref. [la ], but very similar to the OPE behavior 
obtained for the axial double-charmed meson T cc , also described by a scalar-diquark axial-antidiquark 
current [l6j . 

The calculation of the phenomenological side at the hadron level proceeds by writing a dispersion 
relation to the correlation function in Eq. @: 

U phen^ )= f ds P 1 " 1 ™^) + ... (6) 
J s-q z 
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where p phen (s) is the spectral density and the dots represent subtraction terms. The spectral density is 
described, as usual, as a single sharp pole representing the lowest resonance plus a smooth continuum 
representing higher mass states: 



p? h ™( s ) = f 2 S(s-m 2 z )+p cont (s), 
where fz gives the coupling of the current to the meson Z + : 

<0|j|^ + ) = fz- 



(7) 



(8) 



For simplicity, it is assumed that the continuum contribution to the spectral density, p cont (s) in Eq. ([7]), 
vanishes bellow a certain continuum threshold sq. Above this threshold, it is assumed to be given by the 
result obtained with the OPE. Therefore, one uses the ansatz [l7| 



'(•) 



o OPE (s)0(s- So ) 



(9) 



After making a Borel transform to both sides of the sum rule, and transferring the continuum contribu- 
tion to the OPE side, the sum rules for the pseudoscalar meson Z + , up to dimension-eight condensates, 
can be written as: 



fh 



-m%/M 2 



ds e 



-s/M 2 n OPE 



4m - 



where 



m 2 c (qga.Gq)(qq) 
8tt 2 



da exp 



a(l - a)M 2 



1 



a(l - a)M 2 



(10) 



(11) 



To extract the mass mz we take the derivative of Eq. (fTTj|) with respect to 1/M 2 , and divide the result 
by Eq. ([ID]). 

The values used for the quark masses and condensates are [H, [HI]: m c (m c ) = (1.23 ± 0.05) GeV, 
{qq) = -(0.23±0.03) 3 GeV 3 , {qga.Gq) = m^{qq) with m 2 = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 . 

We evaluate the sum rules in the Borel range 2.2 < M 2 < 3.5 GeV 2 , and in the so range 4.8 < ^sq < 5.0 
GeV. 
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FIG. 1: The OPE convergence in the region 2.2 < M 2 < 3.5 GeV 2 for ^ = 4.9 GeV. Perturbative contribution 
(long-dashed line), {g 2 G 2 } contribution (dashed line), {qq} 2 contribution (dotted-line), (qga.Gq){qq) (dot-dashed 
line) and the total contribution (solid line). 



From Fig. [T] we see that we obtain a quite good OPE convergence for M 2 > 2.5 GeV 2 . Therefore, we 
fix the lower value of M 2 in the sum rule window as M^ in = 2.5 GeV 2 . This figure also shows that, 
although there is a change of sign between dimension-six and dimension-eight condensate contributions, 
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as noticed in [19j, both contributions are very small and, therefore, they do not spoil the convergence of 
the OPE. It is also important to mention that the OPE convergence in this case is much better than the 
OPE convergence for the X(3872) meson [lij], and is comparable with the OPE convergence for heavy 
baryons [20l |. 
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FIG. 2: The dashed line shows the relative pole contribution (the pole contribution divided by the total, pole 
plus continuum, contribution) and the solid line shows the relative continuum contribution for y/so = 4.9 GeV. 

To get an upper limit constraint for M 2 we impose that the QCD continuum contribution should 
be smaller than the pole contribution. The comparison between pole and continuum contributions for 
y/so = 4.9 GeV is shown in Fig. [2j From this figure we see that the pole contribution is bigger than the 
continuum for M 2 < 3.05 GeV 2 . The maximum value of M 2 for which this constraint is satisfied depends 
on the value of Sq. The same analysis for the other values of the continuum threshold gives M 2 < 2.85 
GeV 2 for y/so = 4.8 GeV and M 2 < 3.25 GeV 2 for ^/s^ = 5.0 GeV. In our numerical analysis, we shall 
then consider the range of M 2 values from 2.5 GeV 2 until the one allowed by the pole dominance criteria 
given above. 
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FIG. 3: The Z + meson mass as a function of the sum rule parameter (M 2 ) for y/so = 4.9 GeV. The crosses 
indicate the region allowed for the sum rules: the lower limit (cut below 2.5 GeV 2 ) is given by OPE convergence 
requirement and the upper limit by the dominance of the QCD pole contribution. 

In Fig. [31 we show the Z + meson mass, for ^/so = 4.9 GeV, in the relevant sum rule window, with the 
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upper and lower validity limits indicated. From this figure we see that the results are reasonably stable 
as a function of M 2 . 

Using the Borel window, for each value of sq, to evaluate the mass of the Z + meson and then varying 
the value of the continuum threshold in the range y/so = (4.9 ± 0.1) GeV, we arrive at 

(12) 



m z = (4.40 ±0.10) GeV, 



in a very good agreement with the experimental value [4J. 

To check the dependence of our results with the value of the charm quark mass, we fix ^/sq = 4.9 GeV 
and vary the charm quark mass in the range m c = (1.23 ± 0.05) GeV. Using 2.5 < M 2 < 3.05 GeV 2 we 
get: mz — (4.43 ± 0.08) GeV, in agreement with the result in Eq. (fT2|) . Therefore, we conclude that the 
most important sources of uncertainty in our calculation is the value of the continuum threshod and the 
Borel interval. 

We can extend our results to the bottonium analogous state Zbb, considered as a pseudoscalar B* B\ 
molecule, by exchanging the charm quark in Eqs. (H|) to (fTTj) . by the bottom quark. Therefore, in the case 
of the pseudoscalar meson Z bb , using consistently the perturbative MS'-mass mb(mb) — (4.24±0.6) GeV, 
and the continuum threshold in the range 11.2 < y/so < 11.6 GeV, we find a good OPE convergence for 
M 2 > 8.0 GeV 2 . The OPE convergence in this case is even better than the one presented in Fig. 1. We 
also find that the pole contribution is bigger than the continuum contribution for M 2 < 8.25 GeV 2 for 
y/sE < 11.2 GeV, and for M 2 < 9.9 GeV 2 for ^ < 11.6 GeV. For ^ < 11.2 GeV we found no Borel 
window, since M^ ax < 8.0 GeV 2 . 

We find that the results for the Zbb meson mass, in the allowed sum rule window, are very stable as a 
function of M 2 . Taking into account the variation of M 2 and varying so an d to& in the regions indicated 
we get: 



(10.74 ±0.12) GeV 



(13) 



in a very good agreement with the prediction in ref. [2X1 ] . 

For completeness, we also predict the mass of the strange analogous meson Z+ considered as a pseu- 
doscalar D*D\ molecule. The current is obtained by exchanging the d quark in Eq. ([T]) by the s quark. 
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FIG. 4: The OPE convergence for the sum rule for Z s , using yfso = 5.0 GeV. The solid with triangles, long- 
dashed, dashed, solide with circles, dot-dashed and dotted lines give, respectively, m s times the quark condensate, 
four-quark condensate, gluon condensate, m 3 times the mixed condensate, dimension eight condensate and the 
perturbative contributions. The solid line gives the total OPE contribution to the sum rule. 



The expressions obtained in Eqs. ((5|) for (s) and n""^ 9 ^ (g 2 ) should be changed to: 

m 2 c (qq)Q 



(s) 

jjmix(qq) 



47T 2 

mlm,QQ(qq) 
2V 



xfl - Am 2 /s, 



da- 



a(l — a) 



a(l — a)q 2 



1 + 



m 2 — a(l — a)q 2 



(14) 



G 



We get also two new contributions due to the strange quark mass: 



(s) 



3m s 
2V 5 



da 
a 



\— — — -m 2 c (qq) J y((a + (3)m 2 c - aj3s) 

Pmin 



(15) 



Using m s = (0.13 ± 0.03) GeV [221 ]. and the continuum threshold in the range y/so = (5.1 ± 0.1) GeV 
we see, from Fig. 21 that there is a good OPE convergence for M 2 > 2.5 GeV 2 . From Fig. Q] we also see 
that, although proportional to m s , the dimension four condensate m s (qq) (the solid line with triangles) 
is the most important condensate contribution. 

The upper limits for M 2 for each value of are given in Table I. 

Table I: Upper limits in the Borel window for Z 3 . 
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In Fig. [5] we show the relative continumm (solid line) versus pole (dashed line) contribution, for Z s 
using y/so = 5.0 GeV, from where we clearly see that the pole contribution is bigger than the continuum 
contribution for M 2 < 2.80 GeV 2 . 
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FIG. 5: Same as Fig. 2 for Z s using ^ = 5.0 GeV. 



In the case of Z s we get a remarkable stability for the Z s mass, in the allowed sum rule window, as a 
function of M 2 as can be seen by Fig. [6l 



Taking into account the variations on M 



s , m s and m c in the regions indicated above we get: 
(4.70 ±0.06) GeV , (16) 



which is bigger than the D*D\ threshold ~ 4.5 GeV, indicating that this state is probably a very broad 
one and, therefore, it might be very dificult to be seen experimentally. 

In conclucion, we have presented a QCDSR analysis of the two-point functions of the recently observed 
Z+(4430) meson. Due to the closeness of the Z+(4430) mass to the threshold of L>*+(2010)£>i(2420), 
we have followed ref. Q, and have considered the Z + (4430) meson as a D*Di molecule. We have also 
presented a QCDSR study for the analogous mesons and Z s considered as B*Bx molecule and D*Di 
molecule respectively. We find very good OPE convergence for these three four-quark mesons, although 
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FIG. 6: The Z s meson mass as a function of the sum rule parameter (M 2 ) for y/so = 5.0 GeV. The crosses delimit 
the region allowed for the sum rule. 

this is not in general the case for tetraquark states @. We got for Z + a mass in a very good agreement 
with the experimental result. 

In the case of Z s we have obtained a mass bigger than the D*D\ threshold. Therefore, our results 
indicate that the Z s meson is probably very broad. 
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